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DIFFUSION RESTRE1HTE 



6R0CESB DE SSSCOHSTRtrCTION D'UNE ISf&GE X D*TOJ fiEtr 

BE PROJECTIONS "PAR. aPPXiXCATION D'TJNS TRAHTSFOKMEE SIT ' 

ONDCT^ETTE 

DESCRIPTION 
Le sujet de cette invention est un proeadS 
de reconstruction d'une image & partir de jeux de 
projections de cette image et en appliquant una 
transformation en ondelette. 

Les procedes de tomographic consistent a 
examiner un ohj et inanimS ou ua etre vivant par un 
reseau de detecfeedrs mobiles-' qui prend Tine aSrie -de 
vues en toumant auteour de lui. Ces vues sonfc des 
projections de la propariete permettant d'exprimer F 
15 1' image < normal emenfc une absorption d'un rayonnement 

traversant l'objet ou une scintillation d'un corps "Z 
radioactif ingSrS par 1'objet), c«est-S-dire des sotnmes 
de la propriSte le long de lignes traversant 1'objet .ex 
dezinies par la eollimation des dStecteurs, Chaque 
detecteur tnasure un© projection de 1' image a cnaque 
vue - Quand un • norobre suffisant de vues et de 
projections a etS pris, on proeede a 1' inversion des 
resultats pour obfeenir" la valsur de la propriSte a 
cbaque point de l'objet ; cette inversion est 
comparable a 1' inversion d'un systeme d' equations de 
grandes dimensions et pent e'effectuer ouverfcement par 

1 

des methodes algebriques, ou, plus fr€quetmnent r par des 
methodes analytiques par lasquelles on applique des 
operations numSriqu.es successives aux projections, sans 
$0 inverser directemant le systems. Tin grand nombre de 
methodes existent, parmi lesquelles on citejra celle qui 
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PROCEDE DE RECONSTRUCTION D 1 "ONE IMAGE A PARTIR D'UN JEU 
DE PROJECTIONS PAR APPLICATION D'UNE TRANS FORMEE EN 

ONDELETTE 

DESCRIPTION 

Le sujet de cette invention est un procede 
de reconstruction d'une image Sl partir de jeux de 
projections de cette image et en appliquant une 
transformation en ondelette. 

Les procedes de tomographie consistent a 
examiner un objet inanime ou un etre vivant par un 
reseau de detecteurs mobiles qui prend une s^rie de 
vues en tournant autour de lui. Ces vues sont des 
projections de la propriete permettant d'exprimer 
1" image (normalement une absorption d'un rayonnement 
traversant I 1 objet ou une scintillation d'un corps 
radioactif inggre par I'objet), c«est-&-dire des. sommes 
de la propri€te le long de lignes traversant 1 ! objet et 
d<§finies par la collimation des detecteurs. Chaque 
d^tecteur mesure une projection de 1 1 image a chaque 
vue • Quand un nombre suffisant de vues et de 
projections a ete pris, on procdde k 1 1 inversion des 
r£sultats pour obtenir la valeur de la propriety & 
chaque point de 1 1 objet ; cette inversion est 
comparable a l 1 inversion d 1 un systSme d' equations de 
grandes dimensions et peut s ' ef f ectuer ouvertement par 
des methodes alg^briques, ou, plus frequernment , par des 
methodes analytiques par lesquelles on applique des 
operations num£riques successives aux projections sans 
inverser directement le systeme. Un grand nombre de 
mSthodes existent, parmi lesquelles on citera celle qui 
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est decrita dans le brevet fran9ais 2 615 619 qui est 
le premier de cette Squipe de recherche, et le brevet 
francais 2 810 141 plus recent qui pr€sente qraelquea 
res^emblances avec le proc6d§ qui sera dScrit ici* On 
peut au^ei citer un article fdtem»z^U — ^fiarenc©-^ 
cef^lr^tej — rai^de^. grass > Ce3 operations aboutiaeent & ce 
cE * ^ r OIX a PP e H e l a r£fc.roproj action de© rces-uxes, c ! est-&- 
dire an calcul de la valeur de la propriete prise £ 
J^hacwx des points de la ligne de la projection. 

II pent §tre avantageux de travailler avec 
des rSsultats e^riiaea dans le domaine frSqueatiel de 
Fourier , comma le prouve le deuxi&ne document - Des 
trans formSes numerigues d' autre nature ont aussi St£ 
employees. 

On se reportera aux figures 6, 7 et 8 poua? 
una description concrete et schSmatique d'un procSdS de 
prise de meenres ■ toroographiqu.es. Utae source F de 
rayonnement et un isysteme de detection 3 sont mobiles 
le long d'un cafire> axmulaire 2 a des positions 
diaiuetralement cpposees, et un rayonnement conique 
originaire de la source P atteint le systeme de 
detection 3 apres avoir traverse i'ob^et l a etndier. 
La piece essentielle du • systSme de detection 3 est un 
reseau 4 bidimenSionnel de detecteurs 5. Des 
projections R de 1« image tridimensionnelle de l'objet 1 
sont oiesureeg par ceux dee detecteurs s qui sent inclus 
dans un pourtour 15 de 1' "otnbre" de l'objet 1. Un 
nombre important de vues de ce genre est pris a autant 
d'angles Q differents orientation du reseau 4 des 
detecteurs S. II est frequemment considered 
convent ionnel lenient pour les calculs, un reseau fictif 
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est dScrite dans le brevet frangais 2 615 619 qui est 
le premier de cette equipe de recherche, et le brevet 
fran<?ais 2 810 141 plus recent qui pr6sente quelques 
ressemblances avec le proced6 qui sera d€crit ici. On 
peut aussi citer un article (donnez la reference 
complete) ... de Grass. Ces operations aboutissent a ce 
qu'on appelle la rStroproj action des mesures, c'est-a- 
dire au calcul de la valeur de la propriete prise a 
chacun des points de la ligne de la projection. 

II peut Stre avantageux de travailler avec 
des rSsultats exprimes dans le domaine frequent iel de 
Fourier, comme le prouve le deuxieme document. Des 
transform^es numlriques d 1 autre nature ont aussi ete 
. employees . 

On se reportera aux figures 6, 7 et 8 pour 
une description concrete et schematique d' un proc£d6 de 
prise de mesures tomographiques . Une source F de 
rayonnement et un systeme de detection 3 sont mobiles 
le long d'un cadre annulaire 2 a des positions 
diametralement opposees, et un rayonnement conique 
originaire de la source F atteint le systSme de 
detection 3 apres avoir traversS l*objet 1 etudier. 
La pidce essentielle du systeme de detection 3 est un 
rlseau 4 bidimensionnel de detecteurs 5. Des 
projections R de 1 1 image tridimensionnelle de I'objet 1 
sont mesurles par ceux des detecteurs 5 qui sont inclus 
dans un pourtour 15 de l 1 "ombre" de I'objet 1. Un 
nombre important de vues de ce genre est pris a autant 
d 1 angles 9 differents d' orientation du reseau 4 des 
detecteurs 5. II est fr^quemment considSre, 
convent ionnellement pour les calculs, un reseau fictif 
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4 r de dStecteur 5* sur uzx plan de detection Pdet qui 
est parall&le au r£seau 4 rfiel et passe par le centre O 
du cadre .2* Des coordonn6es p et q gout d6f iniee pour^x 
Identifier les detect eurs £ 1 apres leurs lignes et 
5 leurs colonnes. Des calculs de rearrangement , courants 
darvs I'art, permettent de transposer lee mesuree d*un 
r£seau 4 quelconque au reseau fictif 4 1 et d'appliquer 
les aigoritfames de reconstruction & celui-ci r 

Le brevet frangais 2 615 619 sera rappel€ 
10 ici puisqu'il expliquait ett dStail quelles operations 
ntuiieriquee, notairciient celles qui sont dites de filtrage 
et de r$troprojectibn/ permettent d'.obtenir l 1 image de 
I'objet a partir de . ses projections $ des ttiSthodes ~. 
dif f ^rentes exiatenfr fcependaut - 
15 I«e protol&we auguel on s'attaque ici est la 

reduction du teiqps ou du volume des calcule pour r 
inverser le systems des taesures ; * 
On sait gu'il s»agit d'une des limitations ~* 
le© plus graves dans les proc6d€s de tomograpbie, et de " " 
20 nombreuas proeSdSs nouveaux oat St£ congus pour le tuSrae S 
object if que 1 1 invention, dont celui du deuxiSroe brevet 
cit£. 

L'idSe exposSe ici est d' exploiter les 
propri€t€s particuTiSres d 1 tine transfiormee nutnSrique , 

25 dite decomposition e*i oudel^t^a, des P r ^2P^^ ns V olxr 1 
pr£voir des emplacements de pjgo^eoQjc^h nSgligeables ou 
insignifiantes et de ne pas appliquer les ealouls 
d' inversion k ces ^placements , Des perfect ioimements 
off rent encore une plus grande reduction des calculs, 

30 Sous sa forme la E>lue gen&rale, 1 • invention . 

concerne to parocSdS de reconstruction d'une image 
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4 1 de detecteurs 5 1 sur un plan de detection Pdet qui 
est parallele au rgseau 4 reel et passe par le centre O 
du cadre 2. Des coordonnees p et q sont definies pour 
identifier les detecteurs 13 d'aprds leurs lignes et 
leurs colonnes. Des calculs de rearrangement, courants 
dans l'art, permettent de transposer les mesures d'un 
reseau 4 quelconque au r§seau fictif 4 1 et d'appliquer 
les algorithmes de reconstruction a celui-ci. 

Le brevet francpais 2 615 619 sera rappele 
ici puisqu'il expliquait en detail quelles operations 
num§riques, notamment celles qui sont dites de filtrage 
et de retroprojection, permettent d'obtenir 1 1 image de 
I'objet a partir de ses projections ; des methodes 
differentes existent cependant . 

Le probldme auquel on s 1 attaque ici est la 
reduction du temps ou du volume des calculs pour 
inverser le systeme des mesures. 

On sait qu'il s 1 agit d'une des limitations 
les plus graves dans les proced6s de tomographie, et. de 
nombreux procedes nouveaux ont ete congus pour le meme 
object if que 1 1 invention, dont celui du deuxieme brevet 
cite . 

L'idee exposee ici est d 1 exploiter les 
proprietes particulieres d'une transform^e numerique, 
dite decomposition en ondelette, des projections pour 
prevoir des emplacements de projections negligeables ou 
insignif iantes et de ne pas appliquer les calculs 
d 1 inversion 3l ces emplacements. Des perf ectionnements 
off rent encore une plus grande reduction des calculs. 

Sous sa forme la plus generale, 1 1 invention 
concerne un procede de reconstruction d'une image 



4 J*, J^SC J* ^r^L<JC^ 

comprenant ua ggrie cle dScoropositionaC^successives des 
fcisynfi ondelette donnant/des==^agefet^g 
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-de— px o j we t ioxij^broprenaxtfc vune approximafcion St" 
des series Vu.ccelsi.ves de devils nomo1£o<|aes de chacjue 
5 jeu, puis de^^^^^^^»^es_i^^ttes^das jeux de 
projections^* imagettes~~de 1' image a reconstruire et 
une combinaison. dear imagettes a reconstruire par une 
inversion de decomposition, en ondelette donnant ladite 
image, caractfirisS en ce qu'il coniprend/ dans chacune 
10 des series at suceessivement pour des imagettes de 
details ayant ua nombre croissant de pointy une 
• recherche de portions in^ignifiantes , Jugges pour 
depourvues de contenu, et tme recherche de portions 
homothStiques des portions ins ignifi antes pour chacune 
des imagette^t de dStail* qui suivent dang la sSrie ; efc 
eat* ce qu'il _SP^rend ^e-s3iEi^sion_.das rfitropro j ectioi^ss 
pour totifcee^e s porfciona insigniflaiittefct tauten ^1^3 - J - 
r ^ tiono" ^^ S ^ ygy^ 

^ »*S£fc~ antirieur coraprel&i m t e^^l^^lr^ . 

reconstruction d'une iwage par une decomposition en 
ondelette (US - 5 ' 953 388 - A) , qui cependant ^t 
appXiquSe IS, pour reconstruire une portion seulement de 
1* image, en exploitant la propri€t6 de ^localite* de la 
decomposition, q*ii n'est gu&re sensible aux auferes 
portions de 1' image r cpi'cm peut ainsi nggliger dang les 
calculs * 

L 1 invention- sera mainteiiant dScrite plus 
ooacr^tetuent et cc^l&tement en liaison aux figures 
SUivantea : 

30 -la- figure 1 illustre une decomposition 

d r image en ondelette, 
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comprenant un serie de decompositions success ives des 
jeux de projections en ondelette dormant des imagettes 
des jeux de projection, comprenant une approximation et 
des series successives de details homologues de chaque 
jeu, puis de retroprojections des imagettes des jeux de 
projection en imagettes de 1 1 image reconstruire et 
une combinaison des imagettes & reconstruire par une 
inversion de decomposition en ondelette dormant ladite 
image, caract^rise en ce qu'il comprend, dans chacune 
des series et successivement pour des imagettes de 
details ayant un norabre croissant de points, une 
recherche de portions insignif iantes, jugees pour 
depourvues de contenu, et une recherche de portions 
homoth^tiques des portions insignif iantes pour chacune 
des imagettes de details qui suivent dans la serie ; et 
en ce qu'il comprend une omission des retroprojections 
pour toutes les portions insignif iantes et toutes, les 
portions homoth^tiques . 

L'art anterieur comprend un exemple de 
reconstruction d'une image par une decomposition en 
ondelette {US - 5 953 3 88 - A) , qui cependant est 
appliquee la pour reconstruire une portion seulement de 
1 T image, en exploitant la propriety de "localitS" de la 
decomposition, qui n'est guSre sensible aux autres 
portions de 1' image, qu'on peut ainsi n€gliger dans les 
calculs. 

L 1 invention sera maintenant d§crite plus 
concretement et completement en liaison aux figures 
suivantes : 

- la figure 1 illustre une decomposition 
d 1 image en ondelette , 
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- lea figures 2, 3 et 4 illustrent certains 
aspects de 1 • invention, 

- la figure s illustre un organigrarattte 
resunusint une realisation complete du procede, 

- la figure 6 illustre un dispositif de 
principe du proc£d£, 

- et less figures 7 et 8 illustrent plus en 
detail le precede et certaines rotations employees. 

Nous comnenceroos par expos ex la 
transformation d'un signal par une decomposition en 
ondelette. Plueieurs raodeles d'ondelettee existent, qui 
ont en common d'etre comparables a un filtxe passe-fcas. 
Le signal est separe en deux portions dont l'tme, 
associee aux basses frequences/ pent Stre tenue pour 
une approximation- du signal, alors que 1« autre, 
associee aux nautes frequences, exprime plutet see 
details. Xtoe preprints des ondelettes est que les 
portions peuvent coxitenir chacun me moitie dea pointe 
du signal si -bien qu'il n'y" a pas de perte 
d' information par cette decomposition. La decomposition 
peut e.tre faite dans le domaine direct d 1 expression du 
signal ou dans le' domaine de Fourier. 

Dans Te" cas de projections d«un objet 
examine sur un reseafc. generalement bidimensionnel de 
detecteurs, les projections peuvent etre groupies en 
jeux bidimensionnels d'apres deux de leur coordonnees 
(en general p et q sur les axes d'un reseau de 
detecteurs}. Toutefois, corame la figure s 1' illustre, 
on considexe le plus souvent des jeux de projections 
rearrangees sur un reseau fietif de detecteurs. Dans 
l»exemple de la figure 8, on regroupe des projections 
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- les figures 2, 3 et 4 illustrent certains 
aspects de 1 1 invention, 

- la figure 5 illustre un organ igramme 
resutnant une realisation complete du proc6de # 

5 - la figure 6 illustre un dispositif de 

principe du proced6 7 

- et les figures 7 et 8 illustrent plus en 
detail le procedg et certaines rotations employees. 

Nous commencerons par exposer la 
10 transformation d'un signal par une decomposition en 
ondelette. Plusieurs modules d'ondelettes existent, qui 
ont en commun d'etre comparables §. un filtre passe-bas. 
Le signal est separe en deux portions dont l'une, 
associee aux basses frequences, peut etre tenue pour 
15 une approximation du signal, alors que l f autre, 
associee aux hautes frequences, exprime plutot ses 
details. Une propriete des ondelettes est que les 
portions peuvent contenir chacun une moitie des points 
du signal si bien qu'il n'y a pas de perte 
20 d 1 information par cette decomposition. La decomposition 
peut etre faite dans le domaine direct d' expression du 
signal ou dans le domaine de Fourier. 

Dans le cas de projections d ! un objet 
examine sur un reseau general ement bidimensionnel de 
25 detecteurs, les projections peuvent etre groupees en 
jeux bidimensionnels d'apres deux de leur coordonnees 
(en general p et q sur les axes d'un reseau de 
detecteurs), Toutefois, comme la figure 8 l 1 illustre, 
on considere le plus souvent des jeux de projections 
3 0 rearrangges sur un reseau fictif de detecteurs. Dans 
1'exemple de la figure 8, on regroupe des projections 
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TtX originates d'un certain noitfpire de positions Fx 
successives de la source P du ra,yonnement # de maniere 
que les projections aboutissant S, une tneme colonne (S q 
constant) de d£tecteurs 5' du r§seau fictif 4' soient 
5 originaires d'une taeme position Fx, et auesi qua les 
plans de projections soient tous pa.rra.1 Idles jusqu 1 aux 
positions Fx i cm a alors transform^ le probl£me d'une 
geomStrie conique du, rayonneiaent & une gSomStrie en 
Sventails parallSles dont la resolution et plus simple, 
to De plus, le r^seau fictif 4' passe par le centre 0 de 
rotation et appartient done au plan de detection Pdet, 
ce qui faoilite aus#i les calcule* 
^ L f invention pourrait encore &tre appliquSe 

& des rencontructions de coupes de l'objet 1 a l'aide 
d r un r^fieau iRonodimensi.oxinel de dSt^cteurs (tous a 



nsionnex cte aet^cteurs (tons 
m6me coordonnee p4) . On pourrait encore appliquer le 
principe e*pos€ ci-dessus du rearrangement des 
pr ojecta^s^dJ^dne geom£tarie de— ^3?es9S€WKtt& en fiventail 
,plaaa. a^f^oTO^trie .parallSle* ' 
20 fcraitement des laesures rSarrang&es se 

fait en suivant: lee lignes de dgteefceurs 5^d*i yggcgAi 
fictif 4 f , successivamenc pour tous les points. 

Un signal peut; Sere decompose 

successivemeixt en ' ondelette pour donner plusieurs 
25 niveaux de rSsultats* Les nouvelles decompositions ne 
conceraent que la portion de 1' ondelette qui donnait 
1' approximation du signal , la portion ou lee portions 
qui donnaient les details Stanu preservSes. 

Prenons comme exetnple d'objet d'une 
decomposition en ondelette une image formSe de cinq 
cercles dont: un cercle ext£rieur et quatre cercles de 




30 
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Rx originaires d'un certain nombre de positions Fx 
successives de la source F du rayonnement, de maniere 
que les projections about issant a une meme colonne (a q 
constant) de detecteurs 5' du reseau fictif 4 1 soient 
originaires d'une m8me position Fx, et aussi que les 
plans de projections soient tous parall£les jusqu'aux 
positions Fx : on a alors trans forme le problems d f une 
g£ometrie conique du rayonnement & une geom£trie en 
eventails paralleles dont la resolution et plus simple . 
De plus, le reseau fictif 4 1 passe par le centre 0 de 
rotation et appartient done au plan de detection Pdet, 
ce qui facilite aussi les calculs. 

Li 1 invention pourrait encore etre appliqu^e 
a des rencont ructions de coupes de l'objet 1 £ 1 1 aide 
d'un reseau monodimens ionnel de detecteurs (tous a la - 
meme coordonn^e p4) . On pourrait encore appliquer le 
principe expose ci-dessus du rearrangement des 
projections d 1 une geom<§trie en event ail plan a 
g€om€trie parallSle. 

Le traitement des mesures rearrange es se 
fait en suivant les lignes de detecteurs 5 du reseau 
fictif 4 1 , successivement pour tous les points. 

Un signal peut Stre decompose 

successivement en ondelette pour dormer plusieurs 
niveaux de resultats . Les nouvelles decompositions ne 
concernent que la portion de 1' ondelette qui donnait 
1' approximation du signal, la portion ou les portions 
qui donnaient les details etant preservees - 

Prenons comme exempl e d 1 ob j et d ' une 
decomposition en ondelette une image f ormee de cinq 
cercles dont un cercle exterieur et guatre cercles de 



diametres dirferents, tous inscrits. dans le premier. 
L' image da depart comprenait n ac n points , et si on 
applique deux fois, un© decomposition en ondelette de 
cette image selon le prinoipe. precedent, le resultat 
est donne a la figure 1 „ ' 

La decomposition en ondelettes de 1' image 
donne un ensemble d' imagettes, dont trois sont plus 
import antes que les autres, comprenant n/2 x n/2 points 
cnacune, et correspondent auac details horizontaux, aux 
details diagonaux ec aux details verticaux de 1' image 
de depart a granda echelle r elles sont notees 
respectivement Dhl, DdX et Dvl. Les details borizontauac 
de 1> image sont obtenus a partir des projections 
d' angles 6 proobe de o ou jc, les details verticaux a 
partir des projections d' angles 6 proche de ic/ 2 ou 3ic/s/ 
et les details diagonaux a partir des projections 
d» angles intermedia^, res avec les conventions de la 
figure 1. Le reste de 1' image est compose de quatre 
imagettes comprenant n/4 x n/4 points chacune et dont 
trois sont des imagettes de details horizontaux, 
diagonaux et verticaux a plus petite gchelle, notees 
Dn2, Dd2 et Dv2, alors que la derniere imagette est une 
approximation notes de 1' image de depart, 6i on 

appelle <ji et y les fonctions de decomposition d'une 
image ou d'une imagette en ondelette, la fonction . ij> 
donnant 1 » approximation et la fonction y les details, 
les fonctions a appliguer a 1' image de depart pour 
obtenir la decomposition de la figure 1 sont doimees 
par le tableau I. ' 

• : ft (xa) <J> (xa) $ $ 
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diametres diff brents, tous inscrits dans le premier. 
L ' image de depart comprenait n x n points, et si on 
applique deux fois une decomposition en ondelette de 
cette image selon le principe precedent, le rSsultat 
est donne a la figure 1. 

La decomposition en ondelettes de 1' image 
donne un ensemble d' imagettes, dont trois sont plus 
importantes que les autres, comprenant n/2 x n/2 points 
chacune, et correspondent aux details horizontaux, aux 
details diagonaux et aux details verticaux de 1' image 
de depart k grande £chelle ; elles sont not§es 
respectivement Dhl, Ddl et Dvl. Les details horizontaux 
de l 1 image sont obtenus a partir des projections 
d 1 angles 9 proche de o ou K, les details verticaux & 
partir des projections d 1 angles 9 proche de k/ 2 ou 37C/ 2 , 
et les details diagonaux a partir des projections 
d' angles intermediaires avec les conventions de la 
figure 1. Le reste de 1' image est compose de quatre 
imagettes comprenant n/4 x n/4 points chacune et dont 
trois sont des imagettes de details horizontaux, 
diagonaux et verticaux a plus petite echelle, notees 
Dh2, Dd2 et Dv2 , alors que la derniere imagette est une 
approximation not<§e AA de 1' image de depart. Si on 
appelle § et \\f les fonctions de decomposition d'une 
image ou d'une imagette en ondelette, la fonction <|> 
donnant 1 ' approximation et la fonction Y les details, 
les fonctions a appliquer a 1' image de depart pour 
obtenir la decomposition de la figure 1 sont donn^es 
par le tableau I. 

AA : $ * (x 2 ) * (y) 4> 
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2 2 
2 2 

nva , <J> Cxi) ¥ (3^) 

L r invention consiste a effectuer la 
retreprojection sur chacune dee imagettes dee jeux de 
projections decomposers en ondelette,0et a combiner lea 
J5^ettes^gtropro3etees par une inversion de la 
^L^Qd^&oa^. en ©ndelette^ pour obtenir 1' image 
reeherohee. La decomposition en ondelettes est propice 
a diversee simplif ications qui acc61erent beaucoup. les 
calculs. Ces simplifications sont faifces entre la 
IS decomposition et la COmbinaison . 

La premiere d' entre ellee cenceroe des 
filiations qu'on pent Stablir entre des details 
homologues a dee eckelles diffSrentes. On considere 
pour cela des series d» imagettes donnant dee details de 
20 mStne nature. La- • figure 2 (qui illustre une 
decomposition d'une. image reesemblant a celle de la 
figure 1, mais a - trois niveau* de decomposition) 
illustre, pour les trois imagettes des d€tails 
horiaontaux Dbl, Dhfi et Dh3, des portions homologues 
25 Jl, J2 et J3 qui occupent la meme position, et la mSme 
euperficie relative sur chacune de ces imagettes et se 
deduisent ainsi 1'une de 1' autre par une homothetie 
gSometrique dans leurs ■ imagettes . 



Ik 
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Dv2 : * (Xx) $ (x 2 ) (|) (^) Y (^), 
Dd2 , (|> (X!) (J) (x 2 ) ¥ (^) *F (^>, 

Dh2 : (j) (Xx) 4» <x 2 ) ¥ (^) * [&) 

Dvl : (J) (Xi) *¥ (x 2 ) 
Ddl : ¥ (X X ) ¥ (X 2 ) 
Dhl : *F (Xx) * (x 2 ) 

L 1 invention consiste a effectuer la 
r£troprojection sur chacune des imagettes des jeux de 
projections d£composees en ondelette et a combiner les 
imagettes retroprojetees par vine inversion de la 
disposition en ondelette pour obtenir 1' image 
recherchie . La decomposition en ondelettes est propice 
a diverses simplifications qui accelerent beaucoup les 
calculs. Ces simplifications sont faites entre la 
decomposition et la combinaison. 

La premiere d' entre elles concerne des 
filiations qu'on* peut Stablir entre des details 
homologues a des gchelles differentes. On considere 
pour cela des series d 1 imagettes dormant des details de 
mime nature. La figure 2 (qui illustre une 
decomposition d'une image ressemblant a celle de la 
figure 1, mais a trois niveaux de decomposition) 
illustre, pour les trois imagettes des details 
horizontaux Dhl, Dh2 et Dh3 , des portions homologues 
Jl, J2 et J3 qui occupent la m§me position et la meme 
superficie relative sur chacune de ces imagettes et se 
d£duisent ainsi l'une de 1 1 autre par une homoth^tie 
geometrique dans leurs imagettes . 
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On peufc faire I'hypotliSse <iue, pour la 
pl-upart des images rencontrSee eai pratique 
1< exception notainment des images te^turSes) , si vaxe 
portion telle que J3 a un contenu insignif iant, 
5 c'est-a-dire qui ne revele rien par rapport aux valeure 
signif icatives de 1'iroagette, las portions homologues a 
plus grand© echelle, ici J2 et Ji, seront ellea aussi 
insignif iantes - 

Conformement a 1' invention, on. cowmence 

10 done, pour les imagettes des details , 4 retroproQeter 
les details , a l'eofcelle la plus petite, puis les 
details aux ecnelles de plus en plus grandee. Un seuil 
numexigue est applique aux coefficients de l'ondelette, ^ 
Cest-S-dire aux valeurs que prend la transformee dans ^ 

15 l'itnagette considerge,: One valeur inferieure a ce seuil 

donne une portion insignif iante, telle que J3, : & 
Cependant, les portions ineignifiantes des inaagetteg ne ^ 
sont pas reconstruites, c*est-a-dire que les calculs de 
rStropro j ection ne sont pas effectues pour elles. 

20 On procede concrStetnent «n consfcruisant un 

masque nutnSrique avant de retroprojeter l*imagette. Sur 
les details horizontaux Dh, le masque est construct ' 
pour la premiere foig pour l'imagette Dh3. il pr-end une 
valeur egale a o pour les portions insignif iantes 

25 telles que J3 et egale a 1 ailleurs, Les coefficients 
du masque se suiveht dans un crdre determine, par 
exemple ligne apres ligne. Les calculs de 
retroproj ection sont appliques sur 1'imagette 
consideree dans l'ordre des coefficients du masque. 

SO Quand un coefficient est egal a 0, aucun calcul n'est 
effectue pour le point correspondent de l'itnagette, 
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On peut faire I'hypothese que, pour la 
plupart des images rencontr§es en pratique (a 
1' exception notamment des images texturees) f si une 
portion telle que J3 a un contenu insignif iant, 
c'est-a-dire qui ne revele rien par rapport aux valeurs 
signif icatives de 1'imagette, les portions homologues a 
plus grande 6chelle, ici J2 et Jl, seront el les aussi 
insignif iantes . 

Conformement a 1' invention, on commence 
done, pour les imagettes des details, i. retroprojeter 
les details §l l'echelle la plus petite, puis les 
details aux Schelles de plus en plus grande s . Un seuil 
num^rique est applique aux coefficients de l'ondelette, 
e'est-a-dire aux valeurs que prend la transformee dans 
1'imagette consideree. Une valeur inf£rieure a ce seuil 
donne une portion insignif i ante, telle que J3 . 
Cependant, les portions insignif iantes des imagettes ne 
sont pas reconstruites, c'est-&-dire que les calculs de 
retroprojection ne sont pas effectu£s pour elles. 

On procSde concrdtement en construisant un 
masque numerique avant de retroprojeter 1'imagette. Sur 
les details horizontaux Dh, le masque est construit 
pour la premiere fois pour l'imagette Dh3 . II prend une 
valeur egale a 0 pour les portions insignif iantes 
telles que J3 et £gale a 1 ailleurs. Les coefficients 
du masque se suivent dans un ordre determine, par 
exemple ligne apres ligne. Les calculs de 
retroprojection sont appliques sur l'imagette 
consideree dans 1' ordre des coefficients du masque. 
Quand un coefficient est egal a 0, aucun calcul n'est 
effectu€ pour le point correspondant de l'imagette, 
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auquel on attribue une iral^ur nulle jatifst I'itnagette 
retroprojectSe „ 

En abordant 1 ' image suivante des details 
Iiorizontauac (Dh2}, le masque num&rique permet de ne pas 
5 considerer la portion J2, dont les points ne sont pas 
traitSs par l'unit£ de oalcul qui egfiectue la 
rStroprojectipn- Si par exejnple le masque do I'imagette 
Dli3 possSde un coefficient nul & la ligne i et la 
coloxme j f on prSvoit que les quafcxe points des lignee 
10 2i efc 2i+l r et des colonnes 2j et 2j+l de l'ira&gette 
Bh2 auronfc ausei des valour-© insignif iantes - I*es 
calculs de r€troprojection ne iterant pas effectu€s pour 
ees points* 

Ota. masque mimfirique est ainsi consfcruit 

15 pour chaque niveau de dfieoroposition. te masque 
decrxvant I'itnagette Dh2 cotaprendra par convention de© 
coefficients egapsc & 0 pauj? fcoute portion homologue H 
une portion amc coefficients nuls (comma J3) du tnasqiie 
de l'ijnagette correepondante a plus petite Sciielle * 

20 pour determiner les valeurs 1 du 0 des autres points du 
masque de Bh2 f on utilisera encore les cxHqparaisona 
avec le seuil convent ionnel . D' autres portions 
insignifiantes, ayanC des valeurs nulles de masque > 
ponrroAt appara£tre, C'est ce qu'on ^ illustre sur les 

25 imagettes l>vl, Dv2 ' et Dv3, dee d£couipositians 
successive© des details verticauac, L'imagette Dv3 
comprend une portion insignifiante K3 .dont on retrouve 
les homologues K2 et Kl sur les decompositions a plus 
grande ^chell^,, I>'imagatte Dv3 ne coioprend, dans cet 

30 execaple, aucune autre portion insignif iante, mais il a 
St6 possible d'eii trouver trois autres , notSes L2, M2 
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auquel on attribue une valeur nulla dans l'imagette 
retropro j ectee . 

En abordant 1' image suivante des details 
horizontaux (Dh2) , le masque numerique permet de ne pas 
considerer la portion J2, dont les points ne sont pas 
traites par 1' unite de calcul qui - effectue la 
retropro ject ion. Si par exemple le masque de 1 1 imagette 
Dh3 possSde un coefficient nul §. la ligne i et la 
colonne j , on pr£voit que les quatre points des lignes 
2i et 2i+l, et des colonnes 2j et 2j+l de l'imagette 
Dh2 auront aussi des valeurs insignif i antes . Les 
calculs de retroprojection ne seront pas effectues pour 
ces points. 

Un masque numerique est ainsi construit 
pour chaque niveau de decomposition. Le masque 
decrivant l'imagette Dh2 comprendra par convention des 
coefficients egaux a 0 pour toute portion homologue a 
une portion aux coefficients nuls (comme J3) du masque 
de 1 1 imagette correspondante & plus petite echelle ; 
pour determiner les valeurs 1 ou 0 des autres points du 
masque de Dh2 , on utilisera encore les comparaisons 
avec le seuil convent ionnel . D' autres portions 
insignif iantes, ayant des valeurs nulles de masque, 
pourront apparaltre . C'est ce qu'on a illustre sur les 
imagettes Dvl, Dv2 et Dv3, des decompositions 
successives des details verticaux. L'imagette Dv3 
comprend une portion insignif iante K3 dont on retrouve 
les homologues K2 et Kl sur les decompositions a plus 
grande echelle. L'imagette Dv3 ne comprend, dans cet 
exemple, aucune autre portion insignif iante, mais il a 
gte possible d' en trouver trois autres, notees L2 , M2 
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et N2 f sux l'iraagette suivante Pv2. Au niveau de 
decomposition & plus grande Schelle, celui de 
1'imagette Dvi, on ne s'occupera pas retroprojecter 
les points situfis awe portion© LI, Ml at Nl homologues 
de L2, M2 et N2 . 

On va maintenant examiner d'autres 
particularity da procede de 1 1 invention favorables, a 
ime acceleration des calculs- 

La pueioASxe repose sur l'£galite entre la 
transfonaee de Fourier d'une projection de I 'image a un 
angle fixe et la transforraSe de Fourier de 1' image svir 
une droite de mSrae angle passant: par I'origine. 

En ee report ant S la figure 3, ou im jeu de 
projections f"ete converti dans le domaine de Fourier 
pour donaer dee projections de nature feequentielle 
dans le systSme d'aspes notes £i et £2, le support du 
jeu de projection <5Lans l^^d^aine de Fourier coirtprend 
les valeurs comprises entre Wet v 0 pour cowrae pour 
£2. line decomposition en ondelefcte, eomiae celle de la 
. figure l, fait appara£tre 1' approximation encore notSe 
A& dans les frequences les plus basses; autour de 
1'origine £l - £,2 • Q # alors que les details se 
trouvent de part, et d' autre de cette approximation/ & 
des frequences de plus en plus grandes pour les details 
a grande Schelle. 

On a represents une droite passant par 
V origins, en fsisant un angls 0 avec I'axe horizontal 
£l. Cette droite traverse 1' approximation &A, ainsi que 
les imagettes des details verticaux Dvl et Dv2. M^is on 
est assure que les projections faisatxfc cet angle 0 
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et N2, sur 1'imagette suivante Dv2 . Au niveau de 
decomposition a plus grande echelle, celui de 
1'imagette Dvl, on ne s'occupera pas de retropro j ecter 
les points situes aux portions LI, Ml et Nl homologues 
5 de L2, M2 et N2. 

On va maintenant examiner d' autres 
particularity du proced6 de 1 1 invention favorables, a 
tine acceleration des calculs. 

La premiere repose sur l'egalite entre la 
10 transform^ de Fourier d'une projection de 1 1 image a un 
angle fixe et la transform^ de Fourier de 1' image sur 
une droite de m§me angle passant par l'origine. 

En se reportant & la figure 3, oii un jeu de 
projections a et€ converti dans le domaine de Fourier 
15 pour donner des projections de nature f requentielle 
dans le systeme d'axes notes £l et £,2 , le support du 
jeu de projection dans le domaine de Fourier comprend 
les valeurs comprises entre -v 0 et v 0 pour £i comme pour 
£ 2 . Une decomposition en ondelette, comme celle de la 
20 figure 1, fait apparaxtre 1 ' approximation encore notee 
AA dans les frequences les plus basses, autour de 
l'origine £l = C 2 = °/ alors que les details se 
trouvent de part et d' autre de cette approximation, 3. 
des frequences de plus en plus grandes pour les details 
25 a grande echelle. 

On a represente tine droite passant par 

l'origine, en faisant un angle 0 avec 1'axe horizontal 
£l. Cette droite traverse 1 ' approximation AA, ainsi que 
les imagettes des details verticaux Dvl et Dv2 . Mais on 
3 0 est assure que les projections faisant cet angle G 
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seront tout a fait inutiles pour les ealculs de 
rfitroprojection des details diagonaujc Dd et horizontaux 
Dh puisgue la droite d 1 angle 0 passe & l'£cart de leurs 
imagettes, L' application de 1' invention coraprend alora 
5 une selection, pour chacune das categories d' iiuagettes f 
des angles de projections qui seront utilises dans les 
ealculs, Le calcul sur le support des projections est 
glSmentaire* 

ua figure 4 reprend la division d'une image 

XO d€compos6e en ondelette et transpose dans le domaine 
de Fourier. Dne reconstruction 'parfaite peut Stre 
obtenue avec 1 ' ere^>loi d' un nombre dSterminfe de 
projection, d€psndant de la frequence (discretisation) 
de reconstruction de 1* image de l'objet, Ainei, pour 

15 reconstruire 1 4 approximation A&, on peut demofctrer 
qpA'il suffira de c&oisir un nombre de projections, 
panpi celles qui ont €t# faites, correspondent a une 
frequence maafciTCiale VI servant de rayon & un cercle 
circonscrivant la representation fr^quentielle <ie 

20 1 4 approximation dans la decomposition. Be m&ne, les 
details £ petite Schelle seront compl&teraent rendus eg 
utilisant un norabre de projections correspondent & la 
frequence v2 darts le cercle ayan,t oe rayon, dans le 
systSme d'axes £2, et circonscrivant les 

25 representations £r£quentielles du groupe de dStails 
Dv2 f Dd2 et Dh2 a ro§me Schelle* Le support des 
projections dans le domaine de Fourier permet done de 
connaxtre ais&aent les fr&juencee maximal es nScessaires 
pour la retroprojection parfaite des ima^ettes 

3 0 respect ives • 
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seront tout a fait inutiles pour les calculs de 
retroprojection des details diagonaux Dd et horizontaux 
Dh puisque la droite d 1 angle G passe £ 1 1 ecart de leurs 
imagettes . L' application de 1' invention comprend <=ilors 
une selection, pour chacune des categories d' imagettes, 
des angles de projections qui seront utilises dans les 
calculs. Le calcul sur le support des projections est 
eiementaire . 

La figure 4 reprend la division d'une image 
decomposSe en ondelette et transposee dans le domaine 
de Fourier. Une reconstruction parfaite peut etre 
obtenue avec l'emploi d'un nombre determine de 
projections, dependant de la frequence (discretisation) 
de reconstruction de 1 1 image de 1 1 obj et . Ainsi , pour 
reconstruire 1 ' approximation AA, on peut demontrer 
qu'il suffira de choisir un nombre de projections, 
parmi celles qui ont ete faites, correspondant a une 
frequence maxima le vl servant de rayon a un cercle 
circonscrivant la representation f requentielle de 
1 ' approximation AA dans la decomposition. De meme, les 
details a petite echelle seront completement rendus en 
utilisant un nombre de projections correspondant a la 
frequence v2 dans le cercle ayant ce rayon, dans le 
systSme d'axes £l, £2, efc circonscrivant les 
representations f requentielles du groupe de details 
Dv2, Dd2 et Dh2 a meme echelle. Le support des 
projections dans le domaine de Fourier permet done de 
connaitre a i semen t les frequences maximales necessaires 
pour la retroprojection parfaite des imagettes 
respectives . 
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I>a mi^e en ffuvre de l f invention pour 
exploiter ces deux partioularites comprendra done, £ la 
rSferopxoj action de chacune des itnagettes, ucie selection 
des projections utiles pour eette rSfcroprojection, en 
S ecartant les autre a ; et 6ventnellement une restriction 
du nornbre de projections utiles ef fiectiventertt 
exploitSes par le calcul pour n'en garder- que le nombre 
util$. 

La figure 5 est un organigrraxifme de 
10 I'eneeiDble du proced# dScrit ici. 
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La mise en ceuvre de 1 ' invention pour 
exploiter ces deux particularity comprendra done, a la 
retroprojection de chacune des imagettes, une selection 
des projections utiles pour cette retroprojection, en 
ecartant les autres ; et eventuel lenient une restriction 
du nombre de projections utiles ef f ectivement 
exploitees par le calcul pour n'en garder que le nombre 
utile. 

La figure 5 est un organigr amme de 
1 1 ensemble du procedl decrit ici . 
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1. Procede de reconstruction d ! une image a partir de 

de projections de cette ixnagre r comprenant vine 
5 eerie de decompositions success ives des^ jewc de 

0^ 4§HvC A ^tr^J^Z^j^. 

projections en ondelette dotma&t^^ 

^ ^^=fetti&tss ,\ coraprenant des itaage^^ d'une 

approximation (AA) et de£ series successive** {Dd, Dh, 
l>v) de dStaile homologues de chtaque jeu, puis dee 
10 r<6troproj ect iana^g^^i a gS te te e g >deg jeux de projection^ *h 
en imagettes de l 1 image & re coiis t rili re et una 
combiaaison des imagettes a reconstruire par une 
inversion de decomposition dormant ladite image, 
caracterisS en ce- qja'il comprend, dans chacuue des 
IS s&ries et success ivement pour des imagettes de details 
ayanfc un nombre croissant de points , une recherche de 
portions insignif iantes jugees pour depourvues de 
contenu, et une recherche de portions homothetiques des 1 
portions insignif iantea dans chacune des imagettee de 
20 details qui suivenfc dans la serie ; et en ce qu*il 
go^ngrend une omission des rStroprojections pour toutes 
-j, ^^J^Portions insignif iantes efe — i-^hh - 1 gg -^^-^-r^m 

2. Procede de reconstruct ion d'une image suivant la <np^ 
2S revendication 1, earacteriee en ce qu'il comprend une 

selection de regions d' angles (0) des jeuac de 
projections qui sont utilisfis dans les 
retropro j ections , 

3. procSde de reconstruction d'une image suivant la 
30 revendication 2, caracterise en ce qu'il comprend une 
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RE VEND I CAT I ONS 



1. Procede de reconstruction d'une image a partir de 
jeux de projections de cette image,, comprenant une 
serie de decompositions successives des jeux de 
projections en ondelette donnant des imagettes des jeux 
des projections, comprenant des images d ! une 
approximation (AA) et des series successives (Dd, Dh, 
Dv) de details homologues de chaque jeu, puis des 
r<§troprojections des imagettes des jeux de projection 
en imagettes de 1 1 image & reconstruire et une 
combinaison des imagettes & reconstruire par une 
inversion de decomposition donnant ladite image, 
caracterise en ce qu'il cotnprend/ dans chacune des 
series et successivement pour des imagettes de details 
ayant un nombre croissant de points, une recherche de 
portions insignif iante.s jugees pour depourvues de 
contenu, et une recherche de portions homoth£tiques des 
portions insignif iantes dans chacune des imagettes de 
details qui suivent dans la serie ; et en ce qu'il 
comprend une omission des r^troproj ections pour toutes 
les portions insignif iantes et toutes les portions 
homoth^tiques . 

2. Procede de reconstruction d'une image suivant la 
revendication 1, caract£ris£ en ce qu'il comprend une 
selection de regions d 1 angles (6) des jeux de 
projections qui sont utilises dans les 
r^troprojections . 

3. Proced§ de reconstruction d'une image suivant la 
revendication 2, caracterise en ce qu'il comprend une 




selection. d'un norabre de projections qui sont utilisees 
dans la r£troprojection/. 



4. Proc3dS de reconstruction d'une image suivant la 
revendication 2, caraefcerise en ce que la selection des 
regions d 1 angle est rfialisee en fonction d'rm support 
d'une trana£arm6e de Fourier des jeux de projections 

d£cotnpc?egs en ondeldctle. 

5, Proc£d€ de reconstruction d*une image suivant la 
revendioation 3, caracfc&K^e en ce gjie la selection du 
nombre de projections eat rSalisee en f one t ion de 
frequences maxiojales d f un support d'une tra,nsform6e de 
Fourier des jeux de projection^decomposes en ondelette. 
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selection d'un notnbre cle projections qui sont utilisees 
dans la retroprojections . 

4. Procede de reconstruction d'une image suivant la 
revendication 2, caracterise en ce que la selection des 

5 regions d 1 angle est r6alis6e en fonction d ! un support 
d'une transformee de Fourier des jeux de projections 
d^composSs en ondelette. 

5. Procede de reconstruction d'une image suivant la 
revendication 3, caracterise en ce que la selection du 

10 nombre de projections est realis£e en fonction de 
frequences maximales d'un support d'une transform§e de 
Fourier des jeux de projection d§compos6s en ondelette. 
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MESURES (JEUX BIBIMENSIONNELS 
DE PROJECTIONS, OU PREMIERE 
IMAGE DE L'OBJET) 

I 

REAMENAGEMENT DES 
MESURES CONIQUES EN 
MESURES EN EVENTAILS- 
PARALLELES OU DES 
MESURES EN EVENTA1L 
EN MESURES EN PARALLELES 



V 



DECOMPOSITION EN ONDELETTE 
DE L'IMAGE EN IMAGETTES 



SELECTION DES 
PROJECTIONS UTILES 
POUR CHAQUE IMAGETTE 



FIG. 5 



IMAGETTES DES DETAILS A 
LA PLUS PETITE ECHELLE 



DETERMINATION D'UN 
MASQUE DES IMAGETTES 



\/_ 



RETROPROJECTION 
DES IMAGETTES DE 
MEME ECHELLE 



PASSAGE A 
L'ECHELLE 
SUPERIEURE 
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NON 



SYNTHESE DE L'IMAGE 
DE L'OBJET PAR COMBINAISON 
DES IMAGES RETROPROJETEES 
EN ONDELETTE INVERSE 
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